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Abstract 

The azimutal correlations of heavy quarks produced in the high 
energy pp (pp) cohisions are calculated in the framework of QCD with- 
out the usual assumptions of the parton model. The virtual nature of 
the interacting gluons as well as their transverse motion and different 
polarizations are taken into account. We give some predictions for 
the azimutal correlations of charm and beauty hadrons produced at 
Tevatron-collider and LHC. 
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1 INTRODUCTION 



The investigation of the heavy quarks production in the high energy hadron 
processes provides a method for studying the internal structure of hadrons. 
Reahstic estimates of the cross section of the heavy quark production as well 
as their correlations are necessary in order to plan experiments on existing 
and future accelerators. These predictions are obtained usually in the parton 
model framework and depend significantly on the quark and gluon structure 
functions. The last ones are more or less known experimentally from the 
data of HERA, but unknown at the very small values of Bjorken variable 
X < 10"''. However it is just the region that dominates in the heavy quark 
production at the high energies. 

The Gribov-Lipatov-Dokshitzcr-Altarelli-Parisi (GLDAP) evolution equa- 
tion is applied usually to calculate the structure functions. It sums up in the 
leading logarithm approximation (LLA) all the QCD diagram contributions 
proportional to (a^ln 9^)", but it does not take into account the terms pro- 
portional to (q;^ In l/x)". That is why this approximation does not give the 
correct asymptotic behaviour of the structure function in the small x region. 
For the correct description of these phenomena not only the terms of the 
form [as In q^Y l^^ve to be collected in the Feynman diagrams but also the 
terms (a^ In l/x)". 

Another problem that appears at a; ~ is that of the absorption (screen- 
ing) corrections which must stop the increase of the cross section at a; ^ 
in accordance with the unitarity condition. It can be interpreted as the satu- 
ration of parton density. For relatively small virtuality (f < q^{x) the gluon 
structure function behaves as xG{x, q^) ~ (fB?, so the cross section for the 
interaction of the point-like parton with a target, a ~ {l/q^)xG{x,q^) ^ R^, 
obeys the unitarity condition. The quantity qQ{x) can be treated as a new 
infrared-cutoff parameter which plays the role of the typical transverse mo- 
mentum of partons in the parton cascade of the hadron in semihard processes. 
The behaviour of qo{x) was discussed in ref.[6]. It increases with \og{l/x) 
and the values of qQ{x) arc about 2-4 GeV at a; = 0.01 — 0.001. 

The predictions for the cross section of the heavy quark pair production 
are based usually on the parton model calculations [1, 2]. All particles are 
assumed in this model to be on mass shell with the longitudinal component of 
the momentum only (so called coUinear approximation) , and the cross section 
is averaged over two transverse polarizations of the gluons. The virtualities 
(f of the initial partons are taken into account only through their densities. 
The latter are calculated in LLA through GLDAP evolution equation. The 
probabilistic picture of noninteracting partons underlies this way of proceed- 
ing. In the region, where the transverse mass rriT of the produced heavy 
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quark is close to qo{x), the dependence of the amphtude of the subprocess 
99 ~^ QQ) dominating at the high energy, on the virtuahties and polariza- 
tions of the gluons becomes important. The amplitudes of these subprocesses 
should be calculated more accurately than in the parton model. The matrix 
elements of the QCD subprocesses accounting for the virtuahties and po- 
larizations of the gluons are very complicated. We presented them in our 
previous paper [4] for the main and simplest subprocesses, gg — > QQ (~ a1) 
for hadropro duct ion and 7(7 QQ (~ ag) for photo- and electroproduction. 
The contributions of high-order subprocesses can be essential, but our aim 
is to discuss the qualitative difference between our results and the parton 
model predictions. This can be done on the level of the low-order diagrams. 

In this paper we calculate in the same order (~ a^) the azimutal correla- 
tions of the heavy quarks produced in hadron-hadron collisions which are in 
serious disagreement with the conventional parton model predictions (see 
and references therein). 



2 CROSS SECTIONS OF HEAVY FLAVO- 
UR PRODUCTION IN QCD 



The cross section of heavy quark hadroproduction is given schematically by 
the diagrams in Fig.l. The main contribution to the cross section at small x 
is known to come from the gluons. The lower and upper ladder blocks repre- 
sent the two-dimensional gluon distributions ^p{x, qf) and ip{x, q^)-, which are 
the functions of the longitudinal momentum fraction (x and y) of the initial 
hadron and the virtuality of the gluon. Their distribution over x and trans- 
verse momenta in hadron is given in semihard theory |T]] by the function 
ip{x, q^). It differs from the usual function G{x, q^)\ 



Such a definition of (^{x, q^) enables to treat correctly the effects arising from 
the gluon virtuahties. The exact expression for this function can be obtained 
as a solution of the evolution equation which, contrary to the parton model 
case, is nonlinear due to interactions between the partons in the small x 
region. 

In what follows we use Sudakov decomposition for the quark momenta 
Pi,2 through the momenta of the colliding hadrons pa and pb {p\ = — 0) 
and the transverse ones pi^2T'- 




(1) 



Pl,2 = Xi^2PB + yi,2PA + Pl,2T- 



(2) 
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The differential cross sections of heavy quarks hadropro duct ion have the 
form:Q 



dapp 1 1 



- j d^qiTd^q2T5{qiT + q2T - Pit - P2t) 



dyldy^d'^PiTd'^P2T (27r)8 {s] 

X ^^^iq!,yMqlx)\MQQ\'. (3) 

Here s = 2paPb , qi,2T are the gluon transverse momenta and are the 
quark rapidities in the hadron-hadron c.m.s. frame, 



x^=nh^e-y*^, X2=^e-y*2, x = xi + X2 
2/1= ^e^i, y2 = ^eyK 2/ = 2/1 + 2/2 



I Mqq I ^ is the square of the matrix element for the heavy quark pair hadropro- 
duction. 

In LLA kinematics 

Qi ^ ypA + qiT, q2 ^ xpb + q2T, (5) 



so 



QI - -(flT^ (ll - -(l2T- (6) 



2 2 

(The more accurate relations are qf = ~y^' ^2 = ~Ti^' ^lut we are working 
in the kinematics, where x,y ^ 0). 

The matrix element Mqq is calculated in the Born order of QCD without 
standart simplifications of the parton model, since in the small x domain 
there are no grounds for neglecting the transverse momenta of the gluons 
qiT and q2T in comparision with the quark mass and the parameter qo{x). 
In the axial gauge p%A^ = the gluon propagator takes the form D^j,{q) = 
di,u{q)/q^, 

d,M = V - + q''p'B)KPBq). (7) 

For the gluons in t— channel the main contribution comes from the so called 
'nonsense' polarization ^f^^, which can be picked out by decomposing the 
numerator into the longitudinal and the transverse parts: 

5,M) = '^Wa + P>b)Is + 5l ^ 2p%p\/s ^ gl,. (8) 

The other contributions are suppressed by the powers of s. Since the sum of 
the diagrams in Fig.la-lc is gauge invariant in the LLA, the transversality 



"'^We put the argument of as to be equal to the gluon virtuality, which is very close to 
the BLM scheme §1; (see also 0). 
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condition for the ends of gluon line enables one to replace by —q^^jx in 
the expression for g^^. Thus we get 

d,M - -2^, (9) 

X s 

or 

d,M) ^ 2^, (10) 
xys 

if we do such a trick for the vector too. Both these equations for d^J^y can 
be used, but for the form one has to modify slightly the gluon vertex (to 
account for the several ways of gluon emission — see ref. ): 

2 

Te// = [(a^ys - g^y) - gir^sT + 2x iqiTq2T) Pb]- (H) 

■' ■' xys 

As a result the colliding gluons can be treated as aligned ones and their po- 
larization vectors are directed along the transverse momenta. Ultimately, the 
nontrivial azimuthal correlations must arise between the transverse momenta 
Pit and p2T of the heavy quarks. 

Formally, one looses the gauge invariance dealing with off mass shell glu- 
ons (gi, ^2)- Indeed, the new graphs similar to the 'bremsstruhlung' from the 
initial (final) quark line (Fig. Id) may contribute to the central plato rapidity 
region in the covariant Feynman gauge. However this is not the fact. The 
function 0(x, g^) collects with the "semihard" accuracy all terms of the form 
aj(lng^)"(ln(l/x))'" with n + m > k. In this case the triple gluon vertex 
(11) includes effectively all the main (leading logarithmic) contributions of 
the type of Fig. Id 0. For example, the upper part of the graph shown in 
Fig. Id corresponds in terms of the BFKL equation to the t-channel gluon 
reggeization. Therefore the final expression is gauge invariant with the log- 
arithmic accuracy. 

Although the situation considered here seems to be quite opposite to the 
parton model, there is a certain limit in which our formulae can be trans- 
formed into the parton model ones. Consider the pp case and assume that 
the characteristic values of quark momenta pix and p2T are much larger than 
the values of gluon momenta, qiT, (I2T, 

< Pit »< giT > , < P2T »< q2T >, (12) 

and one can keep only the lowest powers of qiT, q2T- It means that we can put 
QiT = '?2T = everywhere in the matrix element Mqq except the vertices. 
Introducing the polar coordinates 

d^iT = Idql^de, (13) 
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(and the same for g2T) and performing angular integration with the help of 
the formula 



r dOrq^^^ql^ = n ql^Sl^y (14) 
Jo 

we obtain 

r deA^^ r ci^2^^M,.M,. = 2vr^|^ |M,„.,p. (15) 
JO y y Jo XX [xyy 

Here Mpart is just the parton model matrix element, since the result is the 
same as that calculated for the real (mass shell) gluons and averaged over 
transverse polarizations. Then we obtain the cross section (^ in the form [2, 
3]: 

da 

dyldy^d'^PiT 

2 1 f as{ql)^{y,qlT)(^s{q2)v{x,qiT) 



l^^-l (i)^y 4V2 7r3 4V2 7r3 ^^^^^^^^ 



Mpartl xG{x,q2)yG{y,q^), (16) 



where s = xys is the mass square of QQ pair. 

On the other hand, the assumption (|12]) is not fulfilled in a more or less 
realistic case. The transverse momenta of produced quarks as well as the 
gluon virtualities (QCD scale values) should be of the order of heavy quark 
masses. 



3 RESULTS OF CALCULATIONS 

Eq.(3) enables to calculate straightforwardly the distributions over the az- 
imutal angle 0, which is defined as the opening angle between the two pro- 
duced heavy quarks projected onto a plane perpendicular to the beam. Let 
us define this plane as xy-plane and let the first heavy quark flights to the x 
direction, i.e. piTx = Pit, ViTy = 0. In this case cos0 = "P^TxIv^t-, and it is 
not a problem to evaluate the distribution over integrating Eq.(3) using, 
say, VEGAS code. 

Since the functions (y9(x, g|) and ^{y^ qf) are unknown at the small values 
of ^2 and qf, we rewrite the integrals in Eq.(3) as 

r dci^5{q.T-piT-P2T)5{y, + y2 - l)^^^iql x)\M\'- = 
Jo qf yo 
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1 IMP 

= AV27r^5ip^T + P2t) Sivi + 1/2 - l)asiQl)xGix, Ql) - (^),,^o (17) 

+ r dq^^5{q,T - Pit - P2t) 5{y, + - l)^^^{ql x)\M\'-, 

where Eq. (^ is usedQ. Thus, we obtain the sum of the three different contri- 
butions: the product of two first terms with i = 1, 2 in the r.h.s. of Eq.(17), 
Wi(</)); the sum of the products of first and second terms, W2{(p)] and, finally, 
the product of second terms, ^3(0). The first contribution, wi(0), is very- 
similar to the conventional LO parton model expression, in which the total 
heavy quark momenta is exactly zero and the angle between quarks is always 
180°. However the angle between two heavy hadrons can be sligtly different 
from this value due to a hadronization process. To take it into account we as- 
sume that in the first contribution, where quarks are produced back-to-back, 
the probability to find the final azimutal angle 4> , 



sm I 



(18) 



pI+Pt 



is determined by the expression 

wiph) = -^e-Pl/Po , (19) 

where po = 0.2 GeV/c and ph is the transverse momentum in the azimu- 
tal plane of the hadron coming from hadronization process. The two last 
contributions, W2{4>) and ^3(0), result in a more or less broad distribution 
over the angle between the produced quarks, so we neglect here their small 
modification in the hadronization. 

It is necessary to repeat that our approach is justified only at small x, 
that is at the high enough initial energies. Unfortunately, the highest energy, 
where the experimental data on the charm azimutal correlations exist, is 
only ^/s = 39 GeV 0. The values of qf and g| in Eq.(3), giving an essential 
contributions to the charm production cross section, are not large enough at 
this energy, so the first contribution in Eq.(17) dominates at the not very 
small Ql value and our predictions coincide practically with the results of 
LO parton model. 

At the higher energies the second and third contributions become large 
enough and lead to some difference with the conventional parton model. The 
predictions of charm pair azimutal correlations at ^/s = 1800 GeV and Ql = 



^Here the value Qq should not be mixed with the function QqIx) discussed in the 
Introduction. 
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4 GeV^ are presented in Fig. 2a. The solid hystogram shows the results for all 
produced charm particles. However in this case some hadronization mech- 
anism can contribute to the azimutal correlations of the charmed hadrons 
with the small relative momenta. To decrease this possible contribution 
the dashed histogram presents the same results for the events, where both 
charmed hadrons have the transverse momenta > PTmin, PTmin = 4 GeV/ c. 

The same results at the energy y/s = 14 TeV are shown in Fig. 2b. Note 
that the predicted azimutal correlations have some energy dependence, they 
become more broad, when the initial energy increases, however this energy 
dependence is weak enough, especially at PTmin = 0. Our (/)-distributions 
seem to be more broad than the results of parton model calculations. 

The similar predictions for the beauty pair azimutal correlations at the 
energies ^/s = 1800 GeV and 14 TeV are presented in Fig. 3 for the two 
values of PTmin, and 8 GeV/c. For both these PTmin values the distributions 
become more broad with the energy increase. 

Contrary to the parton model results, we predict two peaks for all cases, 
the standard one at = 180° and the second peak with the smaller altitude, 
at = 0°, originated from the contribution of the diagram Fig.lc. 

The HERA experimental data on -^2(0;, Q^) at the small x, actually 10~^ < 
X < 10~^, show the singular x-behaviour at the moderate (say, ~ 10^ 
GeV^). Both ZEUS [|1T| as well as HI [|T2| Coll. data can be parametrized as 
with 5 = 0.1^0.25. Such a behaviour at x — is in evident contradiction 



with the unitarity and has to be stopped by a shadow mechanism [ffl, ITB|, P 



To see the possible influence of the shadow effects on the azimutal cor- 
relations we make the simplest assumption that the shadowing modifies the 
gluon distribution in such a way that the real distribution can be written as 

where e -C 1 and xGoix^Q"^) is the bare GRV (HO) gluon distribution [|1^ . 
Both xG{x, Q^) and xGq{x, Q"^) distributions are shown in Fig.4 for the value 
e = 0.01. One can see that the difference between them at the smallest x, 
where the data exist (x ~ 10"''), is about 10 %. 

The calculated results for the azimutal correlations of the heavy flavour 
pairs with "shadowed" gluon distributions, Eq.(20), are presented in Fig.5 
and show that the shadowing does not significantly affect our results. 



4 CONCLUSION 

The above discussion shows that the accounting for the virtual nature of the 
interacting gluons as well as their transverse motion and different polariza- 
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tions results in the qualitative discrepances with the parton model predic- 
tions. Using, say, LO QCD formulae we obtain a considerably more broad 
distribution over the angle between two heavy flavours than the convenient 
NLO parton model (see 0). The reason is that the NLO parton model allows 
the discussed distributions to differ from 5-functions only due to the emission 
of one hard gluon. Our approach, to the contrary, incorporates effectively the 
emission of all evolution gluons via the phenomenological gluon distribution. 

We are grateful to M.G.Ryskin for useful discussions and to E.M. Levin 
who participated at the early stage of the calculations. 

This work is supported in part by INTAS-93-0079, Russian Fund of Fun- 
damental Research (95-2-03145) and the Volkswagen Stifung. 
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Figure captions 



Fig. 1. Low order QCD diagrams for tlie lieavy quark production in 
PP (pp) collisions via gluon-gluon fusion (a-c) and the diagram (d) violating 
formally the gauge invariance, which is restored, however, with logarithmic 
accuracy. 

Fig. 2. The calculated azimutal correlations of the charm pair production 
in pp (pp) collisions at y/s = 1800 GeV (a) and ^/s = 14 TeV (b) for all 
events (solid histograms) and for the events, where both the charm particles 
have pt > Prmin (dashed hystograms). 

Fig. 3. The calculated azimutal correlations of the beauty pair production 
in pp (pp) collisions at ^/s = 1800 GeV (a) and 14 TeV (b) for all events 
(solid histograms) and for the events, where both the beauty particles have 
Pt > PTmin (dashed histograms). 

Fig. 4. Gluon structure function of the nucleon without (solid curve) 
and with (dashed curve) shadow correction, Eq.(19). 

Fig. 5. The calculated azimutal correlations of charm (a) and beauty (b) 
pair production in pp [pp) collisions at a/s = 14 TeV for all events (solid 
histograms) and for the events, where both the heavy flavour particles have 
Pt > PTmin (dashed histograms). 
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